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We address the issue of reducing the resource required to compute information-theoretic quantum 
correlation measures like quantum discord and quantum work deficit in two qubits and higher 
dimensional systems. We show that determination of the quantum correlation measure is possible 
even if we utilize a restricted set of local measurements. We find that the determination allows us to 
obtain a closed form of quantum discord and quantum work deficit for several classes of states, with 
a low error. We show that the computational error caused by the constraint over the complete set of 
local measurements reduces fast with an increase in the size of the restricted set, implying usefulness 
of constrained optimization, especially with the increase of dimensions. We perform quantitative 
analysis to investigate how the error scales with the system size, taking into account a set of plausible 
constructions of the constrained set. Carrying out a comparative study, we show that the resource 
required to optimize quantum work deficit is usually higher than that required for quantum discord. 

We also demonstrate that minimization of quantum discord and quantum work deficit is easier in 
the case of two-qubit mixed states of fixed ranks and with positive partial transpose in comparison 
to the corresponding states having non-positive partial transpose. Applying the methodology to 
quantum spin models, we show that the constrained optimization can be used with advantage in 
analyzing such systems in quantum information-theoretic language. For bound entangled states, we 
show that the error is significantly low when the measurements correspond to the spin observables 
along the three Cartesian coordinates, and thereby we obtain expressions of quantum discord and 
quantum work deficit for these bound entangled states. 


I. INTRODUCTION 

Entanglement [1] as a measure of quantum correla¬ 
tions existing between subsystems of a composite quan¬ 
tum system has been shown to be indispensable in per¬ 
forming several quantum information tasks [2, 3]. To 
deal with challenges such as decoherence due to system- 
environment interaction, entanglement distillation proto¬ 
cols [4] to purify highly entangled states from a collection 
of states with relatively low entanglement have also been 
invented. Parallely, various counter-intuitive findings such 
as substantial non-classical efficiency of quantum states 
with vanishingly small entanglement, and locally indistin¬ 
guishable orthogonal product states [5-8] have motivated 
the search for quantum correlations not belonging to the 
entanglement-separability paradigm. This has led to the 
possibility of introducing more fine-grained quantum cor¬ 
relation measures than entanglement, such as quantum 
discord (QD) [9], quantum work deficit (QWD) [10], and 
various ‘discord-like’ measures [II, 12], opening up a new 
direction of research in quantum information theory. Al¬ 
though establishing a link between the measures of quan¬ 
tum correlations belonging to the two different genres has 
also been tried [13], a decisive result is yet to be found in 
the case of mixed bipartite quantum states. Note, how¬ 
ever, that all these measures reduce to von Neumann en¬ 
tropy of local density matrix for pure states. In recent 
years, the interplay between entanglement distillation and 
quantum correlations such as QD and QWD has been un¬ 
der focus [14]. However, proper understanding of the re¬ 
lation between such measures and distillable as well as 
bound entanglement [15-21] is yet to be achieved. 

There has been a substantial amount of work in deter¬ 
mining QD for various classes of bipartite as well as multi¬ 


partite mixed quantum states [22-24]. A common obser¬ 
vation that stands out from these works is the computa¬ 
tional complexity of the task, due to the optimization over 
a complete set of local measurements involved in its def¬ 
inition [9]. For general quantum states, the optimization 
is often achieved via numerical techniques. It has recently 
been shown that the problem of computing quantum dis¬ 
cord is NP-complete, thereby making the quantity compu¬ 
tationally intractable [25]. The lack of a well-established 
analytic treatment to determine QD has also restricted 
the number of experiments in this topic [26] . Despite con¬ 
siderable efforts to analytically determine QD for general 
two-qubit states [22, 23], a closed form expression exists 
only for the Bell diagonal (BD) states [22]. 

A number of recent numerical studies have shown that 
for a large fraction of a very special class of two-qubit 
states, QD can always be calculated by performing the 
optimization over only a small subset of the complete set 
of local projection measurements [27, 29], thereby reduc¬ 
ing the computational difficulty to a great extent. These 
states are constructed of the three diagonal correlators of 
the correlation matrix, and any one of the three magneti¬ 
zations, which is similar for both the qubits. The subset, 
in the present case, consists of the projection measure¬ 
ments corresponding to the three Pauli matrices, 
and cr’’’. Curiously, the assumption that QD can be opti¬ 
mized over this subset for the entire class of such states 
results only in a small absolute error in the case of those 
states where the assumption is not valid [27, 29]. This 
property also allows one to determine closed form for QD 
for the entire class of such two-qubit states within the 
margin of small absolute error in calculation [27-29]. 

Optimization over such a small subset of the complete 
set of local projection measurements is logical in the situ- 
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ations where a constraint over the allowed local measure¬ 
ments is at work. The knowledge of such a subset may be 
important in quantum estimation theory [30], in relation 
to quantities that are non-linear functions of the quantum 
states, such as the QD, where estimation of the parameter 
values to determine the optimal projection measurement 
depends on the the number of measurements required to 
obtain the desired quantity within a manageable range of 
error. Also, the existence of such subsets has the potential 
to be operationally as well as energetically advantageous 
in experimental determination of the QD and similar mea¬ 
sures for a given quantum state. However, the investiga¬ 
tion of the existence of such special subsets of allowed 
local projectors in the computation of quantum correla¬ 
tions like QD and QWD, as of now, are confined only 
to special classes of quantum states in (8) systems 
[27-29]. A natural question that arises is whether such 
subsets of local projectors can exist for general bipartite 
quantum states. The possible scaling of the absolute er¬ 
ror resulting from the limitation on the number of allowed 
local projectors in the subset is also an interesting issue. 

In this paper, we investigate the issue of simplifying 
the optimization of quantum correlation measures like QD 
and QWD in the case of general two-qubit mixed states of 
different ranks with positive as well as non-positive par¬ 
tial transpose (NPPT) [31] by using a restricted set of 
local projection measurements. We provide a mathemati¬ 
cal description of the optimization of quantum correlation 
measures over a restricted set of allowed local projection 
measurements, and discuss the related statistics of com¬ 
putational error. Using a set of plausible definitions of the 
restricted set, we show that the absolute error, resulting 
due to the constraint over the set of projectors, dies out 
considerably fast. Using the scaling of the error with the 
size of the restricted subset, we demonstrate that even a 
small number of properly chosen projectors can form a re¬ 
stricted set leading to a very small absolute error, thereby 
making the computation of quantum correlation measures 
considerably easier. 

Our method also helps us to find expressions with negli¬ 
gible error for QD and QWD of several important classes 
of quantum states. We demonstrate this in the case of 
two-qubit “A” states [32], which occur, in general, in 
ground or thermal states of several quantum spin mod¬ 
els. Hence, our approach provides a way to study co¬ 
operative phenomena present in such systems with less 
numerical difficulty, as demonstrated here for anisotropic 
XY model in the presence of external transverse field [33] . 
We extend the study to the paradigmatic classes of bound 
entangled (BE) states, where the computation of QD and 
QWD using a special restricted subset is discussed, and 
their analytic forms with small error are determined. The 
results indicate that the error resulting from the restricted 
measurement in the case of states with positive partial 
transpose (PPT) is less compared to the states with non¬ 
positive partial transpose (NPPT). Such constrained op¬ 
timizations can be a powerful tool to study physical quan¬ 
tities in higher dimensions, and to obtain closed forms of 


QD and QWD. 

The paper is organized as follows. In Sec. H, we provide 
a mathematical description of the computation of quan¬ 
tum correlation measures by performing the optimization 
over a constrained set of local projection measurements, 
defining the corresponding absolute error in calculation. 
In Sec. HI, we discuss a set of constructions of the subset 
in relation to the statistics of the error for general two- 
qubit mixed states in the state space. We also consider a 
general two-qubit state in the parameter space, and show 
how symmetry of the state helps in defining the restricted 
subset. We demonstrate how our method can be used 
to determine closed form expressions of quantum corre¬ 
lation measures, and comment on the applicability of the 
method in real physical systems like quantum spin models. 
In particular, we show that the constrained optimization 
technique performs quite well in analyzing the anisotropic 
XY model in transverse field in terms of quantum corre¬ 
lation measures. In Sec. IV, the results on the quantum 
correlations in BE states are presented. Sec. V contains 
the concluding remarks. 

II. DEFINITIONS AND METHODOLOGY 

In this section, after presenting an overview of the quan¬ 
tum correlation measures used in this paper, namely, QD 
and QWD, we introduce the main concepts of the paper, 
i.e., constrained QD as well as QWD by restricting the 
optimization over a small subset of the complete set of 
allowed local measurements. We also discuss the corre¬ 
sponding error generated due to the limitation in mea¬ 
surements. 


A. Quantum discord 

For a bipartite quantum state pas, the QD is defined 
as the minimum difference between two inequivalent defi¬ 
nitions of the quantum mutual information. While one of 
them, given by /(pas) = S {pa)+S {pb)-S (pab), can be 
identified as the “total correlation” of the bipartite quan¬ 
tum system pAs [34], the other definition takes the form 
J_>(pab) = S (pb) — S (pbIpa), which can be argued as a 
measure of classical correlation [9]. Here, pA and ps are 
local density matrices of the subsystems A and B, respec¬ 
tively, S (p) = —Tr [plog 2 p] is the von Neumann entropy 
of the quantum state p, and S{pb\pa) = ThkPkS {pab) 
is the quantum conditional entropy with 

Pab = (llfc 0 7b) PAB (life (g) Is)/pfc, (1) 

and 

Pk = Tr [(Hfe 0 Ib) pab] ■ (2) 

The subscript implies that the measurement, rep¬ 
resented by a complete set of rank-1 projective opera¬ 
tors, {Hfe }, is performed locally on the subsystem A, 
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and /_B is the identity operator defined over the Hilbert 
space of the subsystem B. The QD is thus quantified as 
D = min{/(p^B) — J^{pab)}-, where the minimization is 
performed over the set iSc, the class of all complete sets 
of rank-1 projective operators. One must note here the 
asymmetry embedded in the definition of the QD over the 
interchange of the two subsystems, A and B. Throughout 
this paper, we calculate QD by performing local measure¬ 
ment on the subsystem A. 

B. Quantum work deficit 

Along with the QD, we also consider the QWD [10] 
of a quantum state, defined as the difference between 
the amount of extractable pure states under suitably re¬ 
stricted global and local operations. In the case of a bipar¬ 
tite state pab, the class of global operations, consisting 
of (i) unitary operations, and (ii) dephasing the bipartite 
state by a set of projectors, {H^}, defined on the Hilbert 
space 'H of pab, is called “closed operations” (CO) under 
which the amount of extractable pure states from pab is 
given by Iqo = log 2 dim('H) — S{pab)- And, the class 
of operations consisting of (i) local unitary operations, 
(ii) dephasing by local measurement on the subsystem A, 
and (Hi) communicating the dephased subsystem to the 
other party, B, over a noiseless quantum channel is the 
class of “closed local operations and classical communi¬ 
cation” (CLOCC), under which the extractable amount 
of pure states is /clocc = log 2 dim('H) — mins' {p'ab)- 
Here, p'ab = ThkPkP^B the average quantum state af¬ 
ter the projective measurement {H^} has been performed 
on A, with p\g and pk given by Eqs. (1) and (2), re¬ 
spectively. The minimization in /clocc is achieved over 
Sc- The QWD, W, is given by the difference between the 
quantities Ico (pab) and /clocc (pab)- 

C. Constrained Quantum Correlations: Error in 
Estimation 

We now introduce the physical quantity, which will help 
us to reduce the computational complexity involved in 
evaluation of QD and QWD. In particular, we consider 
the bipartite quantum correlations in the scenario where 
there are restrictions on the complete set of projectors 
defining the local measurement on one of the subsystems. 
Let us assume that constraints on the local measurement 
restrict the class of projection measurements to a subset 
Se {Se ^ Sc), where there are n sets of projection mea¬ 
surements in iSfi. Performing the optimization only over 
the set Se, a “constrained” quantum correlation (CQC), 
Qc, can be defined. We call the subset 5 b as the “ear¬ 
marked” set. Let the actual value of a given quantum 
correlation measure, Q, for a fixed bipartite state, pab, 
be Qa- If the definition of Q involves a minimization, 
log 2 d> Qc> Qa , while a maximization in the definition 
leads to log 2 d > Qa > Qc, where d is the minimum of the 


dimensions among the two parties making up the bipar¬ 
tite state and where we have assumed that log 2 d is the 
maximum value of Qa or Qc- For example, one can define 
the constrained QD (CQD) as 

Dc = m:in[I{pAB) - J^{pab)\, (3) 

•SE 

while the constrained QWD (CQWD) is given by 

Wc = minis( p'ab) - S(pab)]- (4) 

Se 

Note that in general, the earmarked sets for QD and 
QWD, represented by 5^ and S)f respectively, may not 
be identical. 

Evidently, the actual projector for which the quantum 
correlation is optimized may not belong to 5b. Therefore, 
restricting the optimization over the earmarked set gives 
rise to error in estimation of the value of Q for a fixed 
quantum state, pab- Let us denote the absolute error 
occurring due to the optimization over Se, instead of Sc, 
for an arbitrary bipartite state pab, by e, where 

en = \Qc-Qa\, (5) 

with log 2 d > e„ > 0. We call this error as the “volun¬ 
tary” error (VE). One must note that the VE depends on 
the size of the earmarked set, n, as well as the distribution 
of the elements of 5b in the space of projection measure¬ 
ments. When n -A oo, we may (but not necessarily) have 
5b Sc, resulting in Qc -A Qa, whence VE vanishes. 
However, for a finite value of n, we denote the VE by e„. 
Note that £„ also depends on the actual form of the n 
projection measurements in 5b. If e = 0 for a quantum 
state even when the optimization is performed over the 
set 5b, we call the quantum state an “exceptional” state. 

Apart from the quantum information-theoretic mea¬ 
sures having entropic definitions, such as QD and QWD, 
there exists a variety of geometric measures of “discord- 
like” quantum correlations. These measures are based 
on different metrics quantifying the minimum distance of 
the quantum state from the set of all possible classical- 
quantum states [8, 35-40]. Although a collection of dis¬ 
tance metrics have been used to characterize geometric 
measures of quantum correlations, it has been shown that 
out of all the Schatten p-norm distances, only the one- 
norm distance has properties that are rather similar to 
the QD as well as QWD [38-40]. But due to the diffi¬ 
culty in optimizing the measure, analytically closed forms 
of one-norm geometric discord has been obtained only for 
some special types of states, e.g., Bell diagonal states, and 
the X states [40]. Our methodology, along with the tradi¬ 
tional QD and QWD, is applicable also to the geometric 
measures that require an optimization. Motivated by the 
usefulness of QD and QWD in certain quantum proto¬ 
cols [5-7, 41], we choose these measures for the purpose 
of demonstration. 
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III. TWO-QUBIT SYSTEMS 


In the case of a ® system where each of the 
subsystems consists of a single qubit only, the rank-1 pro¬ 
jection measurements are of the form {11^ = U\k){k\U\ 
|fc) = 10), 11)}, where U, a local unitary operator in 5'C/(2), 
can be parametrized using two real parameters, 0, and 
as 


U = 


cos I sin |e*'^ \ 
— sin cos | ) ' 


( 6 ) 


Here, O<0<7r, 0<(/)< 27r, and {|0),|1)} denotes 
the computational basis in C^. Note that 0 and (f) can be 
identified as the azimuthal and the polar angles, respec¬ 
tively, in the Bloch sphere representation of a qubit. Let 
us define a parameter transformation, fg = cos 9, so that 
^ /e < 1- Here and throughout this paper, when¬ 
ever we need to perform an optimization over all rank-1 
projection measurements on a qubit to evaluate a given 
quantum correlation, Q, we choose the parameters fg and 
(j) uniformly in [—1,1] and [0,27r] respectively. An ear¬ 
marked set, in the present case, is equivalent to a subset 
of the complete set of allowed values of fg and cj). 


A. Mixed states with different ranks 


Pri4>) in the case of two-qubit mixed states having NPPT 
or having PPT for both QD and QWD are determined by 
generating 5 x 10® states Haar uniformly for each value of 
r = 2, 3, and 4. We find that in the case of QD as well as 
QWD, both Pf{fg), and P^if)) are uniform distributions 
over the entire ranges of corresponding parameters, fg and 
4>, irrespective of the rank of the state as well as whether 
the state is NPPT or PPT. Note here that for two-qubit 
mixed states of rank-2, almost all states are NPPT while 
the PPT states form a set of measure zero [42], which 
can also be verified numerically. However, in the case of 
r = 3 and 4, non-zero volumes of PPT states are found. 
In 0 systems, all NPPT states are entangled while 
PPT states form the set of separable states [43]. 

The fact that all the rank-1 projection measurements 
are equally probable makes the qualitative features of Qc 
depend only on the geometrical structure of the earmarked 
set, Se, and not on the actual location of the elements of 
the set on the Bloch sphere. In the following, we consider 
four distinct choices of the set Se for both QD and QWD 
in the case of two-qubit mixed states with ranks r = 2,3,4, 
and discuss the corresponding scaling of the average VE. 


Case 1: Se with (fe,4>) distributed over a circle on the 
Bloch sphere 


First, we discuss the case of general two-qubit mixed 
states of different ranks. As described above, an optimal 
set of {fg, (j)) values define the optimal projection measure¬ 
ment for the computation of the fixed quantum correlation 
measure, Q. Let us consider the probability, Pr, that the 
optimal values of the real parameters, fg and </>, for a fixed 
measure of quantum correlation, Q, of a randomly chosen 
two-qubit mixed state of rank r, lie in {fg, fg + dfg), and 
(</),</) -|- dcf), respectively. The fact that the real param¬ 
eters, fg and (j), are independent of each other suggests 
that 

Pr = Pr{fe, fdjdfgdcf = P^{fg)dfgP^{(j))d(j), (7) 

which allows one to investigate the two probability den¬ 
sity functions (PDFs), Pr{fe), and Pr{4’), independently. 
Here, Pr{fg)dfg denotes the probability that irrespective 
of the optimal value of 4>, the optimal value of fg lies be¬ 
tween fg and fg + dfg for the fixed quantum correlation 
measure, Q, calculated for a two-qubit mixed quantum 
state of rank r. A similar definition holds for the proba¬ 
bility Pf{(j))d(l) also. 

In the case of two-qubit systems, an uniform distri¬ 
bution of the projection measurements in the measure¬ 
ment space corresponds to the uniform distribution of the 
{fg,(l)) points on the surface of the Bloch sphere. It is, 
therefore, reasonable to expect that the PDFs, Pf{fg)dfg 
and Pf{4>)d(j), correspond to uniform distributions over the 
allowed ranges of values of fg and (p. To verify this numer¬ 
ically, we consider QD and QWD as the chosen measures 
of quantum correlation. The corresponding P}{fg), and 


We start by constructing the earmarked set with projec¬ 
tion measurements such that the corresponding {fg, p) lies 
on the circle of intersection of a fixed plane with the Bloch 
sphere. Let us assume that the corresponding VE result¬ 
ing from the restricted optimization of Q, for an arbitrary 
two-qubit mixed state of rank r is given by ej), n being the 
size of iS^, and where the corresponding n points on the 
circle are symmetrically placed. Let us also assume that 
the CQC calculated by constrained optimization over the 
set Se with n —?► oo is given by Q'^ and the corresponding 
VE, called the “asymptotic error”, is = \Q'c — Qa\- 
To investigate how fast e” reaches on average with 
increasing n, one must look into the variation of e” — 
against n for different values of r. Here, e” is the average 
value of the VE, e”, and is given by 

f elP:{el)del, (8) 

Jo 

with P”(e”)(ie” being the probability that for an arbitrary 
two-qubit mixed state with rank r, the VE lies between 
e” and e” -|- de” when Qc is calculated over of size n 
defined on the chosen plane. A similar definition holds for 
the average asymptotic VE, e^, and the PDF, 
in the limit n —>• oo, where 

e^o.= f elcPUelr.)d£^o.- (9) 

Jo 

We consider two different ways in which the plane is 
chosen, (a) We fix a value of fg = fg such that the 
corresponding states on the Bloch sphere are given by 
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Figure 1. (Color online.) (a) Schematic representation of the 
earmarked set confined on the circle defined by the intersection 
of the {x,y) plane, fixed by fe = 0, and the Bloch sphere. The 
plane is specified by the eigenbases of and . (b) Schematic 
representation of the earmarked set confined on a set of circles, 
defined by the intersections of a set of planes and the Bloch 
sphere. The planes are considered to be symmetrically placed 
on either side of the (x, y) plane fixed by fe = 0. 

1^) = cos ^ |0) + sin where 4> acts as 

the spanning parameter, (b) In the second option, we fix 
the value oi = (j)' while vary fg with the corresponding 
states 1^). We consider both the scenarios, and investigate 
the scaling of the corresponding average VEs for both QD 
and QWD for arbitrary two-qubit mixed states of different 
ranks. 

(a) Fixed value of fg: Unless otherwise stated, here and 
throughout this paper, we shall fix a plane by assigning a 
value to fg. For the purpose of demonstration, we choose 
fg = 0, fixing the (x, y) plane defined by the eigenbasis of 
the Pauli matrices and , where an arbitrary projec¬ 
tion basis can be written as |^) = -^(|0) -b e"^!!)). The 
earmarked set of size n on the perimeter of the circle of in¬ 
tersection of the (x, y) plane and the Bloch sphere can be 
generated by a set of n projectors of the form U\k){k\U\ 
|fc) = |0), |1), obtained by using a fixed fg = 0, and n eq- 
uispaced divisions of the entire range of (f). Here, the form 
of U is given in Eq. (6). The situation is depicted in Fig. 
1(a). To determine the PDFs, and we 

Haar uniformly generate 5 x 10® random two-qubit mixed 
states of rank r = 2,3, and 4 each. The corresponding 
average VE, ej), and average asymptotic VE, in the 
case of both QD and QWD are determined nsing Eqs. (8) 
and (9) respectively. For both the quantnm correlation 
measnres, the qnantity ej) — is found to have a power- 
law decay with the size, n, of the earmarked set, on the 
log-log scale for all values of r. One can determine the 
functional dependence of £)) over n as 

= + ( 10 ) 

where the fitting constant, k, and the scaling exponent, r, 
are estimated from the numerical data. Fig. 2 shows the 
variations of ej) — with n for both QD and QWD in the 
case of rank-2 two-qubit states. The insets in Fig. 2 shows 
the corresponding variations of ej) — with increasing n 
in the log-log scale. 


QD 

r 

NPPT 

PPT 

3 

K = 1.45 x 10“^ ± 9.97 x 10“^ 

T = 1.93 ± 1.05 x 10-2 
= 9-56 x 10-2 

K = 1.04 x 10-1 ± 7.10 x 10-2 

r = 1.94 ±1.01 x 10-2 
= 6.98 x 10-2 

4 

K = 1.21 x 10-^ ± 8.32 x 10-^ 

T = 1.94 ± 1.02 x 10-2 
= 7.80 x 10-2 

K = 8.67 x 10-2 ± 5.96 x IQ-i' 

T = 1.94 ±9.90 x 10-3 

^4 ^ 5 73 X 10-2 

QWD 

r 

NPPT 

PPT 

3 

K = 1.89 x 10-^ ± 1.30 x 10-2 

T = 1.93 ± 1.13 x 10-2 
= 1.18 x 10-1 

K = 1.52 x 10-1 ± 1.04 x 10-2 

r = 1.93 ±1.08 x 10-2 

^3 ^ 9 g4 X 10-2 

4 

K = 1.50 x 10-1 ± 1.03 x 10-2 

T = 1.93 ± 1.07 x 10-2 
eJU'* = 9.23 x 10-2 

K = 1.20 x 10-1 ± 8.26 x 10-3 

T = 1.94 ± 1.04 x 10-2 
= 7.52 x 10-2 


Table I. Values of the fitting constant, k, scaling exponent, 
r, and asymptotic error, e''^ in the case of NPPT as well as 
PPT two-qubit mixed states of rank r = 3 and 4. Here, QD 
and QWD are considered as quantum correlation measures, 
and the earmarked set is fixed on the circle defined by the 
intersection of a plane fixed by /e = 0, and the Bloch sphere. 
The corresponding power-law variations of with n 

are depicted in Fig. 3. 



5 10 15 20 25 30 35 40 45 50 

n 

Figure 2. (Color online.) Variation of a function 

of n in the case of QD and QWD. (Inset) Linear variation of 
as a function of n in the log-log graph for QD 
and QWD, where the variation is given by Eq. (10). The 
numerical data is represented by points while the fitted curve 
is given by solid lines. The fitting parameters are estimated as 
K = 1.77 X 10”^ ± 1.22 X 10■^ and r = -1.92 ± 1.11 x 10“^ 
with = 1.21 x 10“^ in the case of QD, whereas for QWD, 
K = 2.58 X 10”^ ± 1.78 X 10■^ r = -1.91 ± 1.23 x 10■^ 
and = 1.66 x 10“^. In the main figure, the abscissa is 
dimensionless, while the quantities and are in bits. In 
the inset, the ordinate is in the natural logarithm of 
and the x axis is in the natural logarithm of n. 
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The variations of with n in the log-log scale 

using both QD and QWD with NPPT as well as PPT 
two-qubit mixed states of rank r = 3 and r = 4 are shown 
separately in Fig. 3. The corresponding exponents and 
fitting parameters are quoted in Table I. Note that al¬ 
though the exponent, t, in the case of QD and QWD has 
equal values up to the first decimal point for all the cases 
(for NPPT and PPT states with different ranks), the av¬ 
erage asymptotic VE, is larger in the case of QWD 
in comparison to that for QD. This is reflected in the fact 
that the graph for QWD is above the graph for QD, as 
shown in Fig. 2 and 3. Note also that ej) is less in the 
case of PPT states than the NPPT states, when two-qubit 
states of a fixed rank are considered. 

Note that in the above example, we have fixed the plane 
by fixing fe = 0, which corresponds to a great circle on 
the Bloch sphere. One can also fix a great circle on the 
Bloch sphere by fixing any value of </> in its allowed range 
of values. Earmarked sets defined over any such great cir¬ 
cle on the Bloch sphere have similar scaling properties of 
the average VE as long as the points corresponding to the 
projection measurements belonging to the set Se are dis¬ 
tributed uniformly over the circle. However, for different 
distribution, one can obtain different scaling exponents 
and fitting parameters. We shall shortly discuss one such 
example. 

For fg ^ 0, smaller circles over the Bloch spheres are 
obtained. One can also investigate the scaling behaviour 
of the average VE in the case of earmarked sets having 
elements corresponding to points distributed over these 
smaller circles by using </> as the spanning parameter. 
However, different values for the scaling exponents and 
fitting parameters are obtained as \fg \ —?► 1. 

(b) Fixed value of (p: Next, we consider the scenario 
where the distribution of the {fg, p) points corresponding 
to the projection measurements constituting the set Se on 
the circle of our choice is different (i.e., non-uniform) than 
the previous example. This may happen due to restric¬ 
tions imposed by apparatus during experiment, or other 
relevant physical constraints. As before, we choose the 
plane by fixing <p = 0, thereby confining the earmarked 
set on the great circle representing the intersection of the 
{x,z) plane defined by the eigenbases of and and 
the Bloch sphere. However, to demonstrate the effect of 
such non-uniformity over the scaling parameters, we con¬ 
sider the (/e, (p) points corresponding to the n elements in 
Se by n equal divisions of the entire range of fg. Note that 
the current choice of fg as the spanning parameter leads 
to a different (non-uniform) distribution of the points cor¬ 
responding to the projection measurements in Se on the 
chosen circle. 

Similar to the previous case of fg = 0, one can also 
study the scaling of average VE by defining ej) and cor¬ 
responding to the present case. The variation of ej) — 
with n, as in the previous case, is given by Eq. (10), only 
with different values of fitting constant, scaling exponent, 
and average asymptotic error. In the case of two-qubit 
mixed states of rank-2, the appropriate parameter values 


QD 

r 

NPPT 

PPT 


K = 1.04 x 10-^ ± 5.60 X 10-"‘ 

K = 7.68 X 10-2 ± 4.3 X 10-"^ 


T = 1.47 ± 1.70 x 10-3 

T = 1.48 ± 1.80 X 10-3 

3 

£^3 = 9.54 X 10-2 

£^3 = g 95 X 10-2 


K = 8.95 X 10-2 ± 9.6 X 10-'‘ 

K = 6.23 X 10-2 ± 2.8 X 10-"^ 


T = 1.48 ±3.4 X 10-3 

T = 1.48 ± 1.4 X 10-3 

4 

£^±3 = 7.81 x 10-2 

= 5.79 X 10-2 


Table II. Values of the fitting constant, k, scaling exponent, r, 
and asymptotic error, in the case of NPPT as well as PPT 
two-qubit mixed states of rank r = 3 and 4. Here, quantum 
correlation is quantified by QD, and the earmarked set is fixed 
on the circle of intersection of a plane fixed by <?!> = 0, and the 
Bloch sphere. The n points corresponding to the projection 
measurements in iSb are distributed over the circle by taking 
into account n equispaced division of the range [—1,1] of the 
spanning parameter fg. 


for QD are found to be k = 1.30 x 10 ^ ± 1.08 x 10 
T = 1.47 ± 2.7 X 10-3, with = 1.21 x lO’b In the 
case of states with r = 3 and 4, the values of k, t, and 
in the case of QD are tabulated in Table H. Note that the 
scaling exponents are different from those found in case of 
a fixed fg. However, the fact that the average asymptotic 
VE is less in the case of PPT states compared to that 
of NPPT states remains unchanged even in the present 
scenario. 

One can carry out similar investigation taking QWD 
as the chosen measure of quantum correlation. As in the 
previous case where fg = 0, here also the exponent in the 
case of QWD is found to be same with that for QD up to 
the first decimal place, although the average asymptotic 
error is higher. 


Case 2: Se with (fe,<p) on a collection of circles on the 
Bloch sphere 

We now consider the situation where the projection 
measurements in the earmarked set are such that the cor¬ 
responding {fg,p) are not confined on the perimeter of a 
single fixed disc only, but are lying on the perimeters of a 
set of fixed discs (Fig. 1(b)). As discussed earlier, there 
are several ways in which a disc in the state space can be 
fixed. For the purpose of demonstration, we achieve this 
by fixing the value of fg. The size, n, of the set Se, de¬ 
pends on two quantities: (i) the number, ni, of divisions 
of the allowed range of (p on any one of the discs, and 
(ii) the number, n 2 , of discs that are considered for con¬ 
structing the earmarked set. Evidently, n = nin 2 . Note 
that the value of ni is assumed to be constant for every 
disc, although one may consider, in principle, a varying 
number, nj, such that n = ’^i- 

We demonstrate the situation with an example where 
an arbitrary disc is fixed by fg = fg, and a collection of 
additional discs positioned symmetrically with respect to 
the fixed disc is considered. The fact that the number 
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Figure 3. (Color online.) Variations of ej) — as a function of n, in log-log scale, for r = 3 and 4 in the case of NPPT and 
PPT two-qubit mixed states using QD and QWD. The numerical data is represented by points while the fitted curve is given by 
solid lines. The corresponding values of the fitting constant, k, scaling exponent, r, and asymptotic error, are tabulated in 
Table I. The ordinates of all the figures are in natural logarithm of e)) with and individually being in bits, while the 

abscissa is in natural logarithm of the cardinality of the earmarked set. 


n 2 includes the fixed disc itself implies that n 2 is always 
an odd number. In particular, the discs defining the set 
Se can be marked with different values of fg, given by 
fg = fg±jh, where 0 < j < (n 2 -l)/ 2 , and h = 2 /{n 2 -l). 
The set 5^;, in the present case, approaches the complete 


QD 

QWD 

r 

NPPT 

PPT 

r 

NPPT 

PPT 

2 

ni = Y 

n2 = 9 

- 

2 

721 = y 
712 = 11 

- 

3 

ni = b 

n2 = 8 

ni = 5 

712 = 7 

3 

Til = Y 

712 = 10 

ni = b 

712 = 9 

4 

ni = 5 

n2 = 8 

ni = 5 

712 = 6 

4 

Til = b 

712 = 9 

Til = 5 

712 = 8 


Table III. Values of ni and n 2 that are sufficient to obtain a 
value of average VE of the order of 10~® in the case of NPPT 
and PPT two-qubit mixed states of different ranks, r, for QD 
and QWD. The values correspond to the case where the ear¬ 
marked set is confined on the surface of the Bloch sphere. The 
corresponding depiction is available in Fig. 5. 


set of rank-1 projectors, Sc, when both ni and n 2 tend 
to infinity. Similar to the previous case, one can study 
the variation of average VE with the increase in the size, 
n, of the set Se- However, unlike the previous case, in 
the situation where n —?► oo is consequent to (fg, </>) corre¬ 
sponding to the earmarked set being distributed over the 
entire Bloch sphere, the average asymptotic error, e'^ = 0 
for all quantum correlation measures, since Qc —>■ Qa in 
such cases. 

Fixing ni to be a number for which the average VE 
calculated over an arbitrary disc in the set is considerably 
small, variation of the average VE, ejj, with increasing n 2 , 
where n = nin 2 , can be studied for QD and QWD. Here, 
the average is computed in a similar fashion as given in 
Eq. (8). However, in the present case, is obtained 

by performing optimization of the corresponding quantum 
correlation measure using the current definition of Se- 
We fix fg = 0, and observe that in the case of QD, the 
average VE decreases linearly with increasing n 2 for all 
values of ni (See Fig. 4). This implies that for a fixed 
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ni 


Figure 4. (Color online.) Linear variation of eO as a function 
of n 2 with ni = 10 (n = n\n 2 ) for QD. The earmarked set, 
in this case, is chosen on the perimeters of a set of discs in 
the Bloch sphere, starting from the disc hxed by fg = 0, and 
placing the additional discs symmetrically on either side of 
the fg — 0 disc. The numerical data obtained in each case is 
fitted to a straight line, = mn 2 + c. The fitting parameter 
values are (i) m = —1.10 x 10“®,c = 0.12 (rank-2 states), (ii) 
m = —8.59 X 10“^, c = 9.7 x 10“^ (rank-3 NPPT states), (in) 
m = —6.29 X 10~'*,c = 7.14 x 10“^ (rank-3 PPT states), (iv) 
m = —6.99 X 10~"^, c = 7.96 x 10~^ (rank-4 NPPT states), and 
(v) m = —5.22 X 10“'^,c = 5.90 x 10~^ (rank-4 PPT states). 
The abscissa is dimensionless, while the quantity is in bits. 


ni, the spanning of the surface of the Bloch sphere by the 
perimeters of the discs in the set starting from fg = 0 
occurs linearly with an increase in the value of n 2 . Fig. 4 
depicts the variation of the average VE, ej), as a function 
of 712 for Til = 10, where the linear nature of the variation 
is clearly shown. Similar results hold in the case of QWD 
also. 


Case 3: Se with (fe,cl)) on the surface of Bloch sphere 

Let us now consider the situation where the bases in 
Se are such that the corresponding {fe,4>) are uniformly 
scattered over the entire surface of the Bloch sphere so 
that only rii and 712 equal divisions of the entire range of fg 
and (j), respectively, are allowed, leading to an earmarked 
set of size n = n\n 2 - Similar to the previous case, Se 
S c, and e'^ —i 0 when both ni and 712 approach infinity. 
The average error, in the present case, is determined 
in a similar way as in Eq. (8) with the current description 
of 5^;. From the variation of e)) as a function of ni and 712 
in the case of QD and QWD for two-qubit mixed NPPT 
and PPT states of rank t- = 2, 3 and 4, we estimate the 
minimum number of divisions, 77i, and 772 , required in 
order to converge on a sufficiently low value of e)) (~ 10“^) 
in each case. The corresponding values of ni and 772 are 
given in Table III. It is observed that with an increase in 
the rank of the state, the required size of the earmarked 
set, 77, in order to converge on a sufficiently low value of 
average VE, reduces for both QD and QWD. This feature 
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Figure 5. (Color online.) Variation of e), as a function of 
771 and 772 {n = 711712 ) in the case of NPPT and PPT two- 
qubit mixed states of rank r = 2 and r = 4, where QD and 
QWD are considered as quantum correlation measures. The 
ranges of Tii and 772 marked by A are sufficient to obtain a 
considerably low value of ^ (~ 10~®). The area of the region 
decreases in the case of PPT states compared to NPPT states 
in the case of states with a fixed rank r. Also, the region is 
bigger in the case of QWD compared to that in QD, implying 
a requirement of greater resource in the optimization of QWD. 
The corresponding ranges of 77i and 772 are tabulated in Table 
III. The different shades in the figure correspond to different 
values of e),. All quantities plotted are dimensionless, except 
ej), which is in bits. 


is clearly depicted in Fig. 5 with a reduction in the area of 
the region marked as A. Note also that the required size 
is smaller in the case of PPT states when compared to the 
NPPT states of same rank for a fixed quantum correlation 
measure. 


Case f: Triad 

As the fourth scenario, we consider a very special ear¬ 
marked set, called the “triad”, where the set Se consists 
of only the projection measurements corresponding to the 
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Figure 6 . (Color online.) Profiles of the probability density 
function, P 3 (£ 3 ) against £3 for NPPT and PPT two-qubit 
mixed states of different ranks, r. The distributions are sharply 
peaked around low values of £3 in all cases considered. The 
earmarked set, in this case, is taken to be the triad. The ordi¬ 
nates in all the figures are dimensionless, while the quantities 
£ 3 , r = 2, 3, and 4, are in bits. 


three Pauli operators. This is an extremely restricted ear¬ 
marked set, and one must expect large average VE if Qc 
is calculated for an arbitrary two-qubit state of rank r by 
performing the optimization over the triad. However, in 
Sec. IIIB, we shall show that there exists a large class 
of two-qubit “exceptional” states for which the triad is 
equivalent to Sc, with a vanishing VE. 

Before concluding the discussion on general two-qubit 
mixed states with different ranks, r, we briefly report the 
statistics of the VE, where the subscript n = 3 in 
the present case, denoting the size of the triad. We de¬ 
termine the probability, P^{e'^)d£^, that the value of Qc, 
calculated for a randomly chosen two-qubit mixed state of 
rank r, has the VE between £() and £3 + de'^. To do so, as 
in the previous cases, we Haar uniformly generate 5 x 10® 
NPPT as well as PPT states for each value of r = 2,3, 
and 4. Fig. 6 depicts the variations of normalized 
over the complete range [0,1] of £” for NPPT and PPT 
states of different ranks in the case of both QD as well as 


QWD. It is noteworthy that the distributions are sharply 
peaked in the low-error regions, and for a hxed rank r, the 
probability of finding a PPT state with a very low value 
of £3 is always higher than that for an NPPT state. 


B. Two-qubit states in parameter space 

A general two-qubit state, up to local unitary transfor¬ 
mations [ 22 ], can be written in terms of nine real param¬ 
eters as 

PAB = ^[Ia® Ib + ^ CaaCr^ 0 Crg 

oc—x,y,z 

+ c|Jn0cr^].(ll) 

ct—x,y,z 0—x,y,z 


Here, Caa = (cr“ 0 cr“) are the “classical” correlators given 
by the diagonal elements of the correlation matrix (jcaa | < 
1 ), (Si Ib) and = {Ia 0 cr§) are the single site 

quantities called magnetizations (Ic^ |, jc^j < 1 ), given by 
the elements of the two local Bloch vectors, and I a {Ib) 
is the identity operator on the Hilbert spaces of A {B). 

The maximum rank of the two-qubit state given in Eq. 
(11) can be 4. The probability distribution P[fe,(j)), as 
defined in Eq. (7), is obtained by generating 5 x 10® states 
of the form pab by choosing the diagonal correlators and 
magnetizations randomly from their allowed ranges. Here, 
we drop the subscript r for sake of simplicity. Fig. 7(a) 
(for QD) and (c) (for QWD) show the profiles of P{fg,(j)) 
which have three distinctly high populations around the 
set of values (i) {fg = 0,^ = 0 , 7 r), (ii) {fg = 0,cj) = |), 
and (iii) {fg = ± 1,0 < (j) < tt), which correspond to the 
eigenbasis of cr^, and cr^, respectively, on the Bloch 
sphere. Let us now consider small regions around those 
three high density peaks. They are conveniently marked 
with numbers 1-5 in Figs. 7(b) and (d), and are defined 
as 


1 : -1 < /e < -0.9,0 < (^ < tt, 

2 : 0.9 < /e < 1,0 < (/) < TT, 

^ ■ fg + 4''^ < ,0 < 4> < 

^■fe+{4- 0 < (() < TT, 

5:/|±(</>-|)'<cc^0<<(.<7^, (12) 

with uj = 0.3. In the case of QD, about 56.64% of the total 
number of sample states, given in Eq. 11, are optimized 
in the region marked in Fig. 7(b), while the percentage 
is approximately 42.2% in the case of QWD (Fig. 7(d)). 
Note that these are considerably high fractions, taking 
into account the fact that the area of the marked regions 
combined together is small compared to the entire area of 
the parameter space. 

We also consider the optimization of QD and QWD 
in the case of the two-qubit state given in Eq. (11) by 
confining the earmarked set to a collection of projection 
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Figure 7. (Color online.) The probability distribution landscape, over the plane of (/e,(/I>), in the case of a two-qubit 

state pAB of the form in Eq. (11) in the case of QD (a) and QWD (c). The regions 1-5 are marked on the maps of the 
distribution landscape in the case of QD (b) and QWD (d) so that the corresponding quantum correlation for majority of the 
states is optimized in the marked regions. The definition of the marked regions are given in Eq. (12). The different shades in 
the figure correspond to different values of P[fe,(j))- All quantities plotted are dimensionless, except for (f>, which is in degrees. 


measurements corresponding to a uniformly distributed 
set of points on the entire surface of the Bloch sphere. 
As discussed in Sec. Ill A, we divide the ranges of the 
parameters (f) and fg by ni and n 2 equispaced intervals, 
and perform the minimization of QD and QWD over the 
set Se oi size n = nin 2 - We find that a significantly 
low value of average VE is achieved in the case of QD 
when ni > 2 and n 2 > 4. However, as observed in all the 
previous cases, the required size of the earmarked set, Se, 
to obtain an average VE of same order as in QD, is larger 
(ni > 4,712 > 8) when QWD is taken as the measure of 
quantum correlation, as clearly visible in Fig. 8. 


Special case: Mixed states with fixed magnetizations 

We conclude the section by discussing a special case of 
the two-qubit state given in Eq. (11), where apart from 
the diagonal correlators, any one of the three magneti¬ 


zations is non-zero while the other two magnetizations 
vanish. In particular, we consider a two-qubit state of the 
form 

Pm = Ib + X! c““cr(4 0 erg 

a—x,y,z 

Ib + rUglA® cTg), (13) 

with P = x^y, or 2 . Here, pi) (/3 = x,y,z) is 

the magnetization with pi being the local density matrix 
of the qubit i {i = A,B). The two-qubit X-state [32] of 
the form 


/ ai 0 0 6i \ 

0 02 62 0 I 
0 &2 03 0 I ’ 

\ 61 0 0 04 / 


(14) 


written in the computational basis {jOO), jll), jOl), jlO)} 
is a special case of pm with = z. The matrix elements, 
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Figure 8. (Color online.) Variation of £„ as a function of ni 
and n 2 (n = nin2) in the case of two-qubit state, given in Eq. 
(11), where (a) QD and (b) QWD are considered as quantum 
correlation measures. The ranges of ni and n 2 marked by A 
are sufficient to obtain a considerably low value of ^ (~ 10“^). 
The region is bigger in the case of QWD compared to that in 
QD. The different shades in the figure correspond to different 
values of eO. All quantities plotted are dimensionless, except 
eO being in bits. 


{ui : i = 1, ■ ■ ■ ,4} and {bj : j = 1, 2}, are real numbers, 
and can be considered as functions of the correlators Caa 
{a = x^y,z), and the magnetizations, and m^. The 
importance of two-qubit states of the form given in Eq. 
(14) includes the fact that they are found to occur in the 
quantum information theoretic analyses of several well- 
known quantum spin systems, such as the one-dimensional 
XY model in a transverse field [33], and the XXZ model 
[44]. These models possess certain symmetries that gov¬ 
ern the two-qubit reduced density matrices obtained by 
tracing out all other spins except two chosen spins from 
their ground, thermal, as well as time-evolved states with 
specific time dependence, to have the form given in Eq. 
(14) [12, 45-47]. Since completely analytical forms of QD 
and QWD are not yet available [23] , numerical techniques 
need to be employed, and our methodology show a path 
to handle them analytically. See [28, 29] in this respect. 

For the purpose of demonstration, let us assume that 
the magnetization of pm is along the x direction. Gener¬ 
ating a large number (5 x 10®) of such states by randomly 
choosing the correlators and the magnetizations within 
their allowed ranges of values, and performing extensive 
numerical analysis, we find that for about 99.97% of the 
states, it is enough to perform the optimization over the 
triad {Se consists of the projection measurements corre¬ 
sponding to the three Pauli spin operators) to determine 
actual value of QD. This is due to the fact that these 
99.97% of states are “exceptional”, i.e., the VE, £3 = 0 
for all these states when the optimization is performed 
over Se (as discussed in Sec. 11C). Here, £3 is the VE 
where we drop the superscript r for simplicity. In the 
case of the remaining 0.03% of states, for which the opti¬ 
mal projection measurement does not belong to Se- The 
VE, resulting from the optimization performed over Se, 
is found to be £3 < 2.9088 x 10“®. Similar result has 
been reported in [27], although our numerical findings re¬ 


sult in a different bound. As an example, we consider the 
state Pm defined by = (—1)" x 0.956861, (or c^^) 
= (-1)"* X 0.267575, (orc^'^') = (-l)”+™+i x 0.275867, 

m\ = (-l)P X 0.94976, and = (-1)^ x 0.907559 with 

m, n, and p being integers, even or odd. There exists 
eight such states (corresponding to different values of m, 

n, and p) for which £3 = 2.9088 x 10“®. Amongst the 
set of exceptional states, the optimal projector is for 
about 27.4% states. For the rest of the set of exceptional 
states, QD of the half of them (i.e., 36.3% states) are op¬ 
timized for with projection measurement corresponding 
to cr^ while the rest are optimized for the projector 
corresponding to cr^. 

Interestingly, all of the above results except one remains 
invariant with a change in the direction of magnetization. 
For example, a change in the direction of magnetization 
from X to z results in the optimization of 27.4% of the set 
of exceptional states for projection measurement corre¬ 
sponding to <t\. This highlights an underlying symmetry 
of the state suggesting that the presence of magnetization 
((t“) diminishes the probability of optimization of QD for 
corresponding to cr^, a = x,y,z, provided the state 
Pm belongs to the set of exceptional states. 

Note also that our approach offers a closed form expres¬ 
sion for CQD in the case of two-qubit states of the form 
given in Eq. (14) as 


Dc = S{pa) - S{pae) + min 


S',S'+ 


(15) 


where pA = Tr^jp^^]. The quantities S' and S'^. are 
functions of the matrix elements {ai;i = I,-- - ,4}, and 
{bj',j = 1 , 2 }, and are given by 

S' = (oi -I- 02 ) log 2 (ai -I- 02 ) -b (03 -b 04 ) log 2 (a 3 -b 03 ) 

4 

-^a,log 2 ai, (16) 

i=l 
1 ^ 

-S'i = 1-(17) 

where 

af = 1 + (-I)*y(aa^^a 7 +a 3 ^-a 4 p"+ 4 l 67 ±"&^. 

(18) 

The expression of Dc given in Eq. (15) is exact for the 
exceptional states, and results in a very small absolute 
error in the case of all other two-qubit X states, when 
measurement over qubit A is considered. 

The numerical results for the state pm depends strongly 
on the choice of quantum correlation measures. To 
demonstrate this, we compute QWD instead of QD for 
Pm- Our numerical analysis suggests that irrespective of 
the direction of magnetization, for about 95.51% of the 
two-qubit states of the form pm, the QWD is optimized 
over the triad Se, resulting £3 = 0. These states consti¬ 
tute the set of exceptional states in the case of QWD. For 
the rest 4.49% of states, the assumption that H^p^ £ Se 
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results in a VE, £3 = Wc — Wa < 1.0076 x 10“^, where 
Wc is the CQWD, and Wa is the actual value of QWD. 
Note that the upper bound of the absolute error, in the 
case of the QWD, is much higher than that for QD of pm- 
In contrast to the case of the QD, it is observed that for 
a state Pm with magnetization along, say, the x direction, 
within the set of exceptional states, the QWD for > 53% 
of states - a larger percentage - is optimized for cor¬ 
responding to The rest of the exceptional states, with 
respect to optimization of QWD, are equally distributed 
over the cases where corresponds to (j\ and a\. Sim¬ 
ilar to the case of QD (Eq. (15)), one can obtain a closed 
form expression of the CQWD in the case of two-qubit X 
states as 


Wc = min 


S,S± 


^iPAB)y 


(19) 


where the quantities S and S± are given by S' = 
-ELiOilogaOi, and S± = - log 2 (ai^/ 4 ) 

respectively, with af given in Eq. (18). For the 95.51% 
of states with the form of X states, Wc provides a closed 
form for QWD, when measurement is done on qubit A. 


C. Application: Quantum spin systems 

Now we discuss how the constrained optimization tech¬ 
nique introduced in the paper, and its potential to pro¬ 
vide closed form expresions of quantum correlations with 
small errors, can help in analyzing physical systems. 
For the purpose of demonstration, we choose the spin- 
4 anisotropic quantum XY model in an external trans¬ 
verse field, which has been studied extensively using quan¬ 
tum information theoretic measures [12, 45-47]. Success¬ 
ful laboratory implementation of this model in different 
substrates [48-51] has allowed experimental verification 
of properties of several measures of quantum correlations 
leading to a better understanding of the novel proper¬ 
ties of the model. In this paper, we specifically consider 
variants of the model, viz., (i) the zero-temperature sce¬ 
nario where the model is defined on a lattice of N sites, 
with an external homogeneous transverse field, and (ii) 
the two-qubit representation of the model with staggered 
transverse field, at finite temperature. 


L-qubit system with homogeneous field 

The Hamiltonian describing the anisotropic XY model 
in an external homogeneous transverse field [33] , with pe¬ 
riodic boundary condition, is given by [33, 45] 

H + (1 - 5)o-fcrf+i} Ah'^al, 

( 20 ) 

where J, g (—1 < 5 < 1), and h are the coupling strength, 
the anisotropy parameter, and the strength of the external 


homogeneous transverse magnetic field, respectively. At 
zero temperature and in the thermodynamic limit L —>■ 00 , 
the ground state of the model encounters a quantum phase 
transition (QPT) [52] at A = Ac = 1 (A = J/h) from 
a quantum paramagnetic phase to an antiferomagnetic 
phase [33, 45, 52]. A special case of the model is given by 
the well-known transverse-field Ising model {g = 1). The 
Hamiltonian, H, can be exactly diagonalized by the suc¬ 
cessive applications of the Jordan-Wigner and the Bogoli- 
ubov transformations, and the single-site magnetization, 
mf, and the two-spin correlation functions, c(*“, of the 
spins i and j (i,j € {1, 2, • • • , L}, i ^ j, and a = x,y, z) 
can be determined. Using these parameters, one can ob¬ 
tain the two-spin reduced density matrix, pij, for the 
ground state of the model, which is of the form given 
in Eq. (14), where the matrix elements are functions of 
the single-site magnetization, and two-site spin correlation 
functions. 


For a finite sized spin-chain, the QPT at Ac = 1 for 
|(;| < 1 is detected by a maximum in the variation of 
^ against the tuning parameter. A, which occurs in the 
vicinity of Ac = 1. Here, Q is the measure of quantum cor¬ 
relation, viz., QD and QWD, computed for the nearest- 
neighbour {\i — j\ = 1 ) reduced density matrix ob¬ 
tained from the ground state of the model. With the 
increase in system size, the maximum sharpens and the 
QPT point approaches Ac = 1 as A^ = Ac -f where 

Ac is the value of A at which the maximum of ^ occurs 
for a fixed value of L, a is a dimensionless constant, and 
7 is the scaling index. 


We perform the scaling analysis in the case of the 
transverse-field XY model by using CQD and CQWD as 
observables, and computing their values using Eqs. (15) 
and (19), respectively. Fixing the value of the anisotropy 
parameter at g = 0.5, using CQD, the scaling param¬ 
eters are obtained as a = 0.109 and 7 = 1.215, while 
scaling analysis using the CQWD results in a = 1.031, 
and 7 = 1.515. This indicates a higher value of 7 in the 
case of CQWD, and therefore a better finite-size scaling in 
comparison to CQD. These scaling parameters are consis¬ 
tent with the scaling parameters obtained by performing 
the finite-size scaling analysis using unconstrained opti¬ 
mization for computing QD and QWD numerically. This 
indicates that the CQWD can capture the finite-size scal¬ 
ing features perfectly for the transverse-field XY model 
in the vicinity of the QPT. Therefore, our methodology 
provides a path to explore the quantum cooperative phe¬ 
nomena occuring in quantum spin models, in terms of 
different measures of quantum correlations that involve 
an optimization, in a numerically beneficial way, or ana¬ 
lytically. Fig. 9 provides the log-log plot of the variation 
of \\c — Ac] with L, where Wc (Eq. (19)) is used. 
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Figure 9. (Color online.) Finite size scaling analysis for the 
transverse-field XY model using CQWD, computed from Eq. 
(19), as the observable. The QPT point for a system of size 
L approaches Ac = 1 as = Ac -f aL~'^, where 7 = 1.515 
and a = 1.031. All quantities plotted are dimensionless. The 
abscissa of the figure is in natural logarithm of the number 
of qubits, while the ordinate is in natural logarithm of A, a 
dimensionless quantity. 


Two-qubit system with inhomogeneous transverse field 


We now study the two-qubit anisotropic XY model in 
the presence of staggered transverse magnetic field, rep¬ 
resented by the Hamiltonian 

2 

H2 = J {(1 -b 5 ) 0 -?erf + (1 - + X! 

( 21 ) 


where the strength of the external field on qubit i is 
given by hi {i = 1,2), and all the other symbols have 
their usual meaning. The thermal state of the two- 
qubit system, at a temperature T, is given by pr = 
where the eigenenergies {Ej}, and 
the eigenvectors {|V^j)}, j = 0, • • • ,3 are obtained by di¬ 
agonalizing the Hamiltonian (Eq. (21)), f3 = with 
ks being the Boltzman constant, and P[|a)] = |Q;)(a|. 
Here, Z = tr[e“^^^] is the partition function of the sys¬ 
tem. Similar to the previous example, the form of pr is 
as given in Eq. (14), where the matrix elements are given 
by 


Oi 


02 


0-3 


04 


hi 


1 


1 


cosh(^h+) - h+{hl - 
h^_ cosh(/3h _) + h-{h^_ — 
h‘^_ cosh(/3h_) — h_(h?_ — 


2h?_u L 
1 


1 


2h^u L 


cosh(/3h+) + h_|_(h^ — 


gsinh(/lh+) 

h^u 


62 = — 


sinh(/3h_) 

h_u 


4g2)2 sinh(/3h+) , 
4)2 sinh(/?h_) , 
4)2 sinh(/?h_) , 
4g2)5sinh(/3h+) , 
( 22 ) 



-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 

hj/J 

(a) 


-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 
h./J 
(b) 


Figure 10. (Color online.) Variations of (a) CQD and (b) 
CQWD, as computed from Eqs. (15) and (19), respectively, 
against ^ and with g = 0.5, for the two-qubit system with 
inhomogeneous transverse field at a finite temperature, given 
by J/3 = 1. All quantities plotted are dimensionless, except 
for CQD and CQWD, which are in bits. 


with = 


4g2 


-- 


h- = 


4+ di2^ 


, and 


u = cosh(/3/i+) -b cosh(/?/i_). Using Eqs. (15), (19), and 
(22), one can determine CQD and CQWD for the thermal 
state of the model as functions of the system parameters, 
g, and and the temperature T. 

The variations of CQD and CQWD with ^ and 
as calculated using Eqs. (15) and (19), for g = 0.5 and 
Jj5 = 1, are plotted in Fig. 10. The analytical form in 
Eq. (15) is exact for QD in the present case, while Eq. 
(19) results in a maximum VE, Smax = 6.26 x 10“^ in 
the values of QWD, occurring at the points (^,^) = 
(±1.45, ±0.55), for the given ranges of the system pa¬ 
rameters, viz., |^|,|-^| < 2. However, the qualitative 
features of the variation of QWD with and when 
computed using Eq. (19), are similar to those when the 
QWD is computed via unconstrained optimization. Also, 
if one adopts the constrained optimization technique dis¬ 
cussed in Case 3 of Sec. HI A, e^ax reduces to ~ 10“® 
for ni = 8, and 77,2 = 1. Hence, the value of QWD, with 
negligible error, can be obtained in the case of the thermal 
state of the two-qubit anisotropic XY model in an exter¬ 
nal inhomogeneous magnetic field with very small compu¬ 
tational effort, if the constraints are used appropriately. 
This again proves the usefulness of our methodology. 


IV. RESTRICTED QUANTUM CORRELATIONS 
FOR BOUND ENTANGLED STATES 

From the results discussed in the previous section, it 
comes as a common observation that the VE is less in 
the case of two-qubit PPT states when compared to the 
two-qubit NPPT states of fixed rank for a fixed measure 
of quantum correlation. A natural question arising out 
of the previous discussions is whether the result of low 
VE in the case of PPT states, which until now were all 
separable states, holds even when the state is entangled. 
Since PPT states, if entangled, are always BE, to answer 
this question, one has to look beyond 0 systems 
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and consider quantum states in higher dimensions where 
PPT bound entangled states exist. However, generating 
BE states in higher dimension is itself a non-trivial prob¬ 
lem. Instead, we focus on a number of paradigmatic BE 
states in 0 and 0 systems and investigate the 
properties of the VE in a case-by-case basis. 

It is also observed, from the results reported in the pre¬ 
vious section, that the choice of the triad as the earmarked 
set yields good result in the context of low VE in the case 
of a large fraction of two-qubit states in the parameter 
space. Motivated by this observation, in the following 
calculations, we choose the triad constituted of projection 
measurements corresponding to the three spin-operators, 
{S^, , S^}, in the respective physical system denoted by 
C^, as the earmarked set, Se- Here, S^, (3 = x,y,z, are 
the spin operators for a spin-^^ particle (a system of di¬ 
mension d). Erom now on, in the case of the triad as the 
earmarked set, we discard the subscript ‘3’, and denote 
the VE by e for sake of simplicity. If the dimension is 
2 , the earmarked set is given by the triad constituted of 
the projection measurements corresponding to three Pauli 
matrices. 


A. 0 system 


The first PPT BE state that we consider is in a 
system, and is given by [15] 


Pb = 



/ 

6 

0 

0 

0 

0 

6 

0 

0 

\ 



0 

6 

0 

0 

0 

0 

6 

0 



0 

0 

6 

0 

0 

0 

0 

6 


1 


0 

0 

0 

6 

0 

0 

0 

0 


76+1 


0 

0 

0 

0 

fb 

0 

0 

9b 




6 

0 

0 

0 

0 

6 

0 

0 




0 

6 

0 

0 

0 

0 

6 

0 



V 

0 

0 

6 

0 

9b 

0 

0 

fb 

/ 


(23) 


with 


^ 1 + 6 ^ VT^ 

Jb — r, T 9b — n ! 


(24) 


and 0 < 6 < 1. The state is BE for all allowed values of 
6 . We calculate QD and QWD of the state by performing 
measurement on the qubit. In Eig. 11(a), we plot the 
CQD, Dc, and the actual QD, Dq, as functions of the state 
parameter, b. At 6 = 0.15, QD shows a sudden change in 
its variation. The inset of Eig. 11(a) shows the variation 
of the relative error, e, with 6, exhibiting a discontinuous 
jump from zero to a non-zero value at b = 0.15. Note 
that for 6 < 0.15, £r = 0, indicating B^j G Se, which, in 
this case, corresponds to For 6 > 0.15, the maximum 
value of e, although non-zero, is of the order of 10“^, and 
monotonically decreases to zero at 6 = 1. Therefore, a 
restriction of the minimization of QD over the triad, in 
this region, is advantageous. Eor b > 0.15, we observe that 
minimization in obtaining the value of Dc is achieved from 
the projection measurement corresponding to <t^ . Using 




Figure 11. (Color online.) Variation of QD (a) and QWD 
(b) as a function of b in the case of the PPT BE state pb in 
0 C'* systems. (Inset) Variation of the corresponding VE, 
e, as a function of b. The error jumps from zero to a non-zero 
value at 6 = 0.15 in the case of QD, while at 6 = 0.22 in the 
case of QWD. The quantities Dc, Da, Wc, Wa, and e are in 
bits, while b is dimensionless. 


this information, one can obtain a closed form expression 
for Dc as 


Dc = S{pb) - S{pb) + min 


Sl,S2 


(25) 


where S{p^) is the von Neumann entropy of the local 
density matrix of the qubit part. The quantities Si and 
S 2 are functions of the state parameter, b, given by 


Si = 


1 + 761 
2 


1 + 96 + (1 + 36) log2(l + 36) — 261og2 6 


log2 0 


2=1 

2 


= l0g2 Dj + T-j l0g2 T-j), 


(26) 

(27) 


io'=l 
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where the quantities Ci, r^, and r' are given by 


0 = l + 6 +(-l)Vl- 6 ^ 

^ 4(1 1 7b) + (~ 1 ) Vl - + (-iF^z), 

^ 4(1 + 7b) + (-l)Vl - + (-l)^a;'), 

(28) 

with = 2 [l — 36 + 126^ + (—1)*(1 — 36)-\/l — 6 ^J, and 

uj'^ = 2 [l + 6 + 86 ^ + (—1)®(1 + 6 )Vl — 6 ^]. The expres¬ 
sion in Eq. (25) provides the actual value of QD upto an 
absolute error e ^ 10 “^. 

The variations of QWD and CQWD, as functions of 6 , 
are depicted in Fig. 11(b) with the inset demonstrating 
the corresponding variations of e. Similar to the case of 
QD, the optimal measurement observable is for 6 < 
0 . 22 , and for 6 > 0 . 22 , which can be used to determine 
analytic expression of Wc as 

Wc = min[;5i,5'2]-I-S'(/9b), (29) 


where the functions Si and S 2 are given by 
1 


= 


l-h76L 
, 2 


1 -I- 6 -I- (1 -I- 76) log 2 (l -I- 76) - 661og2 6 


“9 0 


(30) 


2=1 

2 


S 2 = -t: X! log 2 Tij + Ty log 2 rb - - rb). 


*. 1=1 


(31) 


The quantities Q, r^-, and rb are given in Eq. (28). In the 
case of QWD, the point at which the sudden change takes 
place (6 = 0.22) is different from that for QD (6 = 0.15). 
Note that the behaviors of both constrained as well as 
actual QWD as functions of 6 are similar to the behaviors 
of the respective varieties of QD, and the VE, in both 
cases, also show similar variations. The maximum value 
of e, in the case of QWD, is also of the order of 10“^. 




Figure 12. (Color online.) Variation of QD (a) and QWD 
(b) against a in the case of the PPT BE state pa in C® ® C® 
systems. (Inset) Variation of the corresponding VE, e, as a 
function of a. The error is maximum at a = 0.23 in the case 
of QD, and at o = 0.32 in the case of QWD. The quantities 
Dc, Da, Wc, Wa, and e are in bits, while a is dimensionless. 


B. C® 0 C® systems 


Case 1 


As an example of a BE state in a 0 systems, we 
first consider the following state [15]: 


Pa = 


1 

8 n -\-1 


/aOOOaO 0 0 a \ 

OaOOOOOOO 
OOaOOOOOO 
OOOaOO 0 0 0 

aOOOaO 0 0 a 

OOOOOaOOO 
000000 /' 0 5' 

OOOOOOOaO 
V a 0 0 0 a 0 5 ' 0 /' / 


(32) 


where 0 < a < 1 , and the functions /' = fb=a and 
g'a = gb=a, with fb and gb given by Eq. (24). Note 
that similar to the state, given in Eq. (23), the state 
is BE for the entire range of a. However, unlike the 
(0 state discussed in Sec. IV A, the local measure¬ 
ment must be performed on a subsystem. We consider 
the set {jo), | 1 ), | 2 )} as the computational basis for each 
subsystem of dimension d = 3. The triad, in this case, 
consists of the measurements corresponding to the observ¬ 
ables S^, S^, and S^, where S^,S = x,y,z, are the spin 
operators for a spin-1 particle. Fig. 12(a) demonstrates 
the variations of QD and CQD as functions of a for the 
entire range [0,1]. Note that the VE (shown in the inset) 
is non-zero for the entire range of a, attaining its maxi¬ 
mum at a = 0.23. The optimization of Dc, is obtained 
for the measurement corresponding to the observable 
when a < 0.23, while for a > 0.23, the optimal measure- 
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(a) 


Figure 13. (Color online.) Variation of Dc and Da for pa, 
as functions of the state parameter a, when measurement is 
performed over the subsystem B. (Inset) Variation of VE, e, 
as a function of a. The quantities Da, Da, and e are in bits, 
while a is dimensionless. 


ment observable is . Similar results are obtained in the 
case of QWD and CQWD (Fig. 12(b)) where the max¬ 
imum VE occurs at a = 0.32. In both the cases (QD 
and QWD), the maximum VE is of the order of ~ 10“^, 
which is much higher compared to the same in previous 
examples. 

An interesting observation comes from swapping the 
subsystem over which the local measurement is performed 
since the state pa is asymmetric over an exchange of the 
subsystems A and B. If one optimizes QD of pa over 
a complete set of local measurements performed on the 
party B instead of A, it is observed that the VE reduces 
drastically in comparison to the same obtained when mea¬ 
surement is performed over the party A. The variation of 
the corresponding Be and Da along with VE against the 
state parameter a is given in Fig. 13. The VE attains 
a non-zero value (from the zero value) at a = 0 . 6 , while 
a sudden change in the variation profile is observed at 
a = 0.665. This change is due to a transition of the opti¬ 
mal measurement observable from to aX a = 0.665. 
Note that the maximum value of VE in the region a > 0.6 
is ~ 10“^, in contrast to the previous case. Note also 
that analytical expressions for CQD and CQWD can be 
obtained in a procedure similar as in the previous case. 


Case 2 


Let us take another example of a PPT BE state in 0 
C^, given by [17]: 


2 a 5 — a 

Pa = + y£'-s H- —P- 


(33) 



Figure 14. (Color online.) Variation of QD (a) and the corre¬ 
sponding VE (b) as functions of a in the case of the PPT BE 
state Qa, given in Eq. (33). The quantities Da, Da, and e are 
in bits, while a is dimensionless. 


where p+ = (|01)(01| -b |12)(12| -b |20)(20|)/3, = 

(|10)(10| + |21)(21| + |02)(02|)/3, ^ E?=o 1**)^ and 

0 < a < 5. The state is separable for 2 < a < 3, BE for 
3 < a < 4 and 1 < a < 2, whi[g)distillable for 4 < a < 5 
and 0 < a < 1 [17]. Fig. 14(a) depicts the variation 
of QD and CQD as functions of the state parameter a 
over its entire range. Note that the QD as well as the 
CQD remains constant over the range of a in which en¬ 
tanglement is distillable, whereas both of them attains a 
minimum at a = 2.5 in the separable region. In the BE 
region, 1 < a < 2, the value of QD as well as CQD re¬ 
mains constant up to a = 1.36, and then decreases for 
increasing a in the region 1.36 < a < 2. Besides, in the 
BE region, 3 < a < 4, the value of QD as well as CQD 
increases with increasing a up to a = 3.64, and then be¬ 
comes constant. The corresponding e is plotted against 
a in Fig. 14(b). Note that the maximum VE is com¬ 
mitted in the separable region (of the order of ^ 10“^), 
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whereas it is smaller (of the order of ^ 10“^) in the BE 
region. Clearly, the maximum VE is relatively higher in 
the present case in comparison to the ® BE state, 
Pb (Eq. (23)), but considerably lower when compared to 
the 0 BE state, Pa (Eq. (32)). When QD is con¬ 
stant, the optimal measurement observable is while it 
changes to either or when value of QD increases or 
decreases. Note that in the latter case, choice of or 
as the optimal measurement observable are equivalent in 
the context of optimizing QD. Analytic expression for Dc 
can be obtained using the above analysis, as shown in the 
previous cases. 

We conclude the discussion on the state Qa by pointing 
out that the QWD of the state coincides with the QD 
since Qa has maximally mixed marginals, i.e., the local 
density operator of each of the parties is proportional to 
identity in [10]. 

V. CONCLUDING REMARKS 

To summarize, we have addressed the question whether 
the computational complexity of information theoretic 
measures of quantum correlations such as quantum dis¬ 
cord and quantum work deficit can be reduced by per¬ 
forming the optimization involved over a constrained sub¬ 
set of local projectors instead of the complete set. We have 
considered four plausible constructions of such a restricted 
set, and shown that the average absolute error, in the case 
of two-qubit mixed states with different ranks, dies down 
fast with the increase in the size of the set. Quantitative 
investigation of the reduction of error with the increase 
in the size of the restricted set has been performed with 
a comparative study between quantum discord and quan¬ 
tum work deficit, and the corresponding scaling exponents 
have been estimated. We have also considered a general 


two-qubit state up to local unitary transformation, and 
have shown that the computation of measures like quan¬ 
tum discord and quantum work deficit can be made con¬ 
siderably easier by carefully choosing the restricted set of 
projectors. We have also pointed out that insight about 
constructing the constrained set can be gathered from the 
probability distribution of the optimizing parameters. 

If a very special restricted set consisting of the projec¬ 
tion measurements corresponding to only the three Pauli 
matrices is considered, we have demonstrated that the 
state space of a two-qubit system contains a large frac¬ 
tion of states for which exact minimizations of quantum 
discord and quantum work deficit are obtained only on 
this set, resulting vanishing error. We have also pointed 
out that this feature can be utilized to obtain closed-form 
expressions of quantum correlation measures up to small 
error for some special classes of states, which can be used 
to study physical systems such as quantum spin mod¬ 
els. The usefulness of this methodology has been demon¬ 
strated in the finite size scaling analysis of the well-known 
transverse-field XY model at zero temperature, and for 
the thermal state of a two-qubit XY model in an external 
inhomogeneous transverse field. Moreover, we have found 
that the absolute error in the value of quantum correla¬ 
tion calculated using the constrained set in the case of 
two-qubit PPT as well as PPT BE states is low compared 
to the NPPT states. The investigations show that these 
measures can be obtained with high accuracy even when 
restrictions in the optimizations involved in their defini¬ 
tions are employed, thereby reducing the computational 
difficulties in their evaluation. Although we have inves¬ 
tigated only information theoretic measures, this study 
gives rise to a possibility to overcome challenges in the 
computation of other quantum correlation measures that 
involve optimization over some set. 
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